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Abstract 

We present an approximate analytical solution of the Cahn-Hilliard equation 
describing the coalescence during a first order phase transition. We have iden- 
tified all the intermediate profiles, stationary solutions of the noiseless Cahn- 
Hilliard equation. Using properties of the soliton lattices, periodic solutions of 
the Ginzburg-Landau equation, we have construct a family of ansatz describ- 
ing continuously the process of destabilization and period doubling predicted 
in Langer's self similar scenario [1]. 
Pacs numbers : 05.45.Yv, 47.20. Ky, 47.54. +r 

I. INTRODUCTION 

When a homogenous system departs suddenly from equilibrium, the fluctuations around 
the initial ground state are linearly amplified and the homogenous phase can for example 
spontaneously separate into two different more stable states. The interfaces which delimit 
the numerous resulting monophasic domains will interact with each other, either giving rise 
to formation of a complex pattern, or merging into a single interface when the domains of the 
same state slowly coalesce, minimizing the overall interfacial energy. It results then in only 
two well separated domains. This process of first order phase transition arises particularly 
for binary mixtures [2] alloys [3], or vapor condensation [4]. 

It can either initiate via a nucleation process, where the homogenous state is put sud- 
denly in a metastable configuration, and an energy barrier has to be crossed before the 
transition appears. Or via a spinodal decomposition when the system is led into a linearly 
unstable configuration. In this latter case, the leading instability selects a modulation of the 
order parameter at a well defined length scale, which will grows and, due to non-linearities, 
saturates. The resulting pattern is composed of well defined interfaces delimiting domains 
containing one of the two stable phases. Remarkably, this intermediate stage conserves the 
modulation width, and the resulting stationary pattern is of almost the same length scale as 
the one selected initially [5,6]. The dynamics finally ends with a much slower, self- inhibiting 
process, dominated by the interactions between the interfaces. The different regions of each 
phase coalesce in the so-called Ostwald ripening where the number of domains diminishes 
whereas their typical size increases. The asymptotic state is decomposed into two domains. 



1 



one for each phase. In this article, by spinodal decomposition, we refer to the first stage of 
the dynamics only, while coarsening will denote the second stage. Although this coarsening 
dynamics is in fact already present, its influence can be often neglected during the first stage 
of the process. 

Hillert [3], Cahn and Hilliard [7] have proposed a model equation describing the segre- 
gation for a binary mixture. This model, known as the Cahn-Hilliard equation (C-H later 
on), belongs to the Model B class in Hohenberg and Halperin's classification [8]. It is a 
standard model for phase transition with conserved quantities and has applications to phase 
transition in liquid crystals [9], segregation of granular mixtures in a rotating drum [10] , or 
formation of sand ripples [11,12]. It is a partial differential equation to which a conservative 
noise is added to account for thermal fluctuations [13]. 

Figure (1) shows snapshots of a numerical integration of the (C-H) dynamics which rep- 
resents the full phase transition process after a quench in temperature. Thermal fluctuations 

were present in the initial conditions, but have been omitted in the dynamics. The three 
main stages of the spinodal decomposition described above are clearly distinguished: first, 
from Fig. (1 (a)) to Fig. (1 (b)), we observe the selection of a typical length scale for the 
modulations, then, from Fig. (1 (b)) to Fig. (1 (c)), the non-linear growth and its satu- 
ration. We note that the number of peaks has been almost conserved between these two 
configurations. On the contrary, during the coarsening dynamics observed between Fig. (1 
(c)) and Fig. (1 (d)) the typical length of the pattern is increasing while on the other hand 
the amplitude of the modulation slowly growth to reach its asymptotic value. 

An important activity has been devoted to the description of the dynamics of phase 
transition, using both statistical methods or numerical simulations (for a review see [14]). 
The late stage of the spinodal decomposition where the coarsening dynamics dominates 
exhibits "dynamical scaling" : the dynamics presents a self-similar evolution where time 
enters only through a length scale L{t), associated with a typical length of the domains 
or the rate of decay of the inhomogeneities. For instance, scaling arguments and stability 
criteria give the law L(t) ~ t^^^ for spatial dimensions greater than one and a logarithmic 
behavior for one dimension in the case of the (C-H) equation [14]. 

This last stage, as observed in two-dimensional demixion of copolymer [15] and as sug- 
gested initially by Langcr [1], can be described as a self similar process of synchronous fusion 
and evaporation of domains. This observation motivated our work and the aim of this arti- 
cle is to present a one dimensional ansatz describing continuously the coalescence process. 
This ansatz is in the form of a one parameter family of symmetric profiles which interpo- 
lates between two stationary states composed of homogeneous domains of length A/2 and 
A. It allows to realize a self similar sequence of coalescence process in ID, starting from the 
regular micro phase separated states issued from the non-linear saturation of the spinodal 
decomposition dynamics and ending with the single interface which characterize the infinite 
time, thermodynamic stable state. 

The paper is organized as followed: first, we present a brief review on general properties 
of phase segregations and on the (C-H) model, mainly to fix the notations. We will repro- 
duce briefly the original derivation by Cahn and Hilliard, restricting ourselves to the one 
dimensional case. In part III, we present a family of symmetric solutions of the Ginzburg- 
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Landau equation which is used to stTidy the dynamics of spinodal decomposition and to 
determine all the symmetric stationary state of the (C-H) dynamics. Then in part IV, the 
main original part of this work, we introduce a non-symmetric family of solutions of the 
(G-L) equation which is used to construct a continuous interpolation between two consecu- 
tive symmetric stationary states. After a study of the energy landscape associated with this 
ansatz, we finally discuss the numerical accuracy of our calculations. In the conclusion, we 
justify the hypothesis we have made and compare the suggested scenario with coalescence 
in real systems. 



II. THE CAHN-HILLIARD MODEL 

The Cahn-Hilliard theory is a modified diffusion equation; it is a continuous conservative 
model for the scalar order parameter which reads in its dimensionless form: 

^ (r, t) = V2(|$ + 2$3 - V^$) + e (r, t) . (1) 

The real order parameter can correspond to the dimensionless magnetization in Ising fer- 
romagnet, to the fluctuation of density of a fluid around its mean value during a phase 
separation or to the concentration in some region around r of one of the components of a 
binary solution, e is the dimensionless control parameter of the system ; it is often identifled 
to the reduced temperature (e = where Tc is the critical temperature of the phase 

-t c 

transition). This equation, flrst derived by Cahn and Hilliard [7], has also been retrieved by 

Langer [1] from microscopic considerations. As written, the (C-H) equation does account 
for thermal fluctuations present in the system through a random white noise ^ (r, i), whose 
amplitude is proportional to the square root of the temperature of the system. 

The homogeneous stationary solutions for the noiseless (C-H) equation are extrema of the 
effective Ginzburg-Landau potential V{^) = -|- $^ (G-L later on). For positive s, there 
is only one homogenous solution $ = which is linearly stable; for negative e, the stationary 
solution $ = undergoes a pitchfork bifurcation and three stationary solutions exist. $ = 
is still a stationary solution, but it is now linearly unstable ; two other symmetric solutions 
$ = i"^^ ^-re stable and have the same free energy F — —e'^/32. Thus, a first order 
transition can be experienced by quenching the system suddenly from a positive reduced 
temperature e to a negative one. Spinodal decomposition is the resulting dynamics. 

The stability of the solution $ = can be studied by linearizing equation (1) around 
$ = (i.e. neglecting the non linear term $^); considering $ as a sum of Fourier modes: 

$(r, t)^Yl </'qe^'*-^+'^* (2) 
q 

where (pg is the Fourier coefficient at t — 0, we obtain for the amplification factor a{q) : 

<^) = -{<f + \W (3) 

It shows immediately that $ = is linearly stable for £ > while a band of Fourier modes 
are unstable for negative £, since a{q} > for < g < sj {—e/2). Moreover, the most 
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unstable mode is for qc-H = a/— £/2(with (Tmax = fg)- This wave number of maximum 
amplification factor will dominate the first stage of the dynamics; in particular, it explains 
why the modulations appear at length scales close to Xc-h = '^t^/Qc-Hi the associated wave 
length. Later on, interfaces separating each domain interact through coalescence dynamics, 
causing < A > to change slowly toward higher values [17,5] 

We will now use known results on non-homogeneous solutions of the (G-L) equation to 
study both the saturation of the spinodal decomposition and the coalescence. 



III. STATIONARY STATES OF THE CAHN-HILLIARD DYNAMICS 



A. Symmetric Soliton Lattice Solutions 



For £ < 0, there exists a stationary solution of the one dimensional (C-H) that relies 
the two homogenous phases $ = i"^^ 



$(x) = ^tanh(^). (4) 

Such a monotonic solution describes a continuum interface between the two stable homo- 
geneous phases, and corresponds to the thermodynamically stable solution that ends the 
phase transition dynamics. But this is a particular member of a one parameter family of 
stationary solutions of the (G-L) equation 

+ 2*=^ - V^* = (5) 

These solutions, the so-called sohton-lattice solutions, are : 



= A;ASn(|, k) with ^ = A'^ = ^ 2^^^! (6) 

where Sn(a;, k) is the Jacobian elliptic function sine-amplitude, or cnoidal mode. This family 
of solutions is parametrized by e* and by the modulus k e [0, 1], or "segregation parameter" . 
These solutions describe periodic patterns of periods 

A = AK{k)i, where K{k) = [ ' , (7) 

Jo VI - A;2 sin^ t 

is the complete Jacobian eUiptic integral of the first kind. Together with k, it characterizes 
the segregation, defined as the ratio between the size of the homogeneous domains, 0.5 x A, 
and the width of the interface separating them, 2 x ^. The equation (7) and the relation 
^ = A~^, enable to rewrite this family as : 

^ , , 4K(k)-k^ AKik) 

= A^ ^' ^ ^ 

This family of profiles (or alternating interfaces) can be obtain exactly as a periodic sum of 
single solitons and antisolitons [18] 

y (-1)" tanh(7rs(x - n)) = ^M^l^Sn(x, k) with s = and k'^ = I - k'' (9) 

^-^ ns K[k') 
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B. Antsatz for the Spinodal Decomposition Dynamics 



The preceding family of profiles can be used to explore the spinodal decomposition 
dynamics. It can be associated with a micro phase separation, locally limited by the finite 
diffusion coefficient. For k — 1, Sn(x, 1) = tanh(x), we recover the usual single interface 
solution (4), of width 2/^/\e\; it is associated with a one soliton solution and corresponds 
to a strong, or macroscopic segregation. Note that K{1) diverges ; the solution 



*i,.(2:) = ^tanh(^x). (10) 

is thus the limit of infinite s, when the solitons, entering in relation (9), are far apart one each 
others. In the opposite limit (weak segregation regime), it describes a sinusoidal modulation 

limfe_o*fc,£(x) = k^l^sm{^j^x) = /c^ sin(^a;) = kqsm{qx) (11) 

It will correspond to the Fourier mode g = ^ of the initial white noise, with an arbitrary 
small amplitude u = kq. Since experiences, numerical simulations and linear stability anal- 
ysis show that A, the spatial period of the pattern is constant during the whole spinodal 
decomposition process, we choose A to coincide with the most instable wave length obtain 
with the Cahn Hilliard linear approach, A = Xc-h = J^^^o ^ where Eq is the quench tem- 
perature. Thus, we obtain a one parameter family of profiles '$*{x,k) — ^fe,Ac-if(^) which 
describe very well both the linear growth and the saturation . The dynamics is now reduced 
to the time evolution of the single free parameter : k(t). Using equations (6) and (7), we 
find that A, k and e are related to one another through the state equation 

sik)^-2il + k^)(^^y. (12) 

So, this implicit equation tells us that if we fixe A = Xc-h, the dynamics can also be reduced 
to the evolution of e(t). 

Given a periodic function $ (obtained either from experimental data or numerical sim- 
ulation of equation (5)) at time t, the ansatz assumes that it corresponds to a soliton 
lattice of the same period: i.e., there exists k{t) such that 'ifk,\c-Hi^) ^ ^{x,t) for each 
time t. For this purpose, we have developed three different algorithms, taking advan- 
tage of the general properties of the family of solutions \&* {x, k) : cither, k can be de- 
duced both from the amplitude of the oscillation equals to 4kK{k) / A, or from the relation 
k =1 — ((\E'(A/2, /c)/\I'(A/4, A;))^ — 1)^; thirdly, a straightforward computation relates k to 
the ratio of the two first Fourier modes of the soliton lattice 'if*{x,k). We have observed 
that the three methods show in general similar results within an error of one percent. 

In this approach, e{t) can be then interpreted as a fictitious temperature or "local tem- 
perature" of the domains: it is the temperature extracted from the profile at a given time, 
using the correspondence between e and k of equation (12). For instance, at t = 0, the 
amplitude is small and we find that k{0) = — > and thus £*{0) = Stt^/A^, different a 
priori from Eq (£*(0) = ^ for A = Xc-h)- 
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Somehow, the dynamics of (C-H) can be projected at first order onto a dynamics along the 
sub-family k) = '^k,\c-n{-'^')^ which can be considered as an attractor of the solutions, 

i.e. the density profile of the system will evolve with time, staying always close to a function 
^*(a;, k). And using a solubility condition, it is possible to compute the full non linear part 
of this dynamics, the saturation of the spinodal decomposition, which leads the system in a 
well defined stationary state [6] . 



C. Saturations of the Spinodal Decomposition Dynamics 

According to the previous interpretation of the parameter e, as e{t = 0) = the system 
is initially out of equilibrium. The dynamics will saturate when this fictitious temperature 
will reach the real thermodynamic one, i.e. the quench temperature So] that is, using 
equation of state (12) for A = Xc-h, when k — solution of the implicit equation : 

2(1 + k'^)K{k'f ^ ^ that is F = 0.687 (13) 

Note that in this case, the width of the interface, which was initially, just after the quench, 

proportional to ^^j^ has now become proportional to "y=^^(^ = — ^^z^ - the 

segregation has slightly increased. Using linear stability analysis, Langer has shown that 
the profile thus obtained, '^*{x,k'^) = \E'(a;, fc*, Ac-if), is destroyed by stochastic thermal 
fluctuations and he has identified the most instable mode as an " antiferro" mode, leading to 
a period doubling. The result of this destabilization is another profile of alternate interface, 
where the length of the domains is now :A = 2\c-h = This means that the new 

stationary profile is given by ^(a;, k{, 2Xc-h), where kf is solution of the implicit relation 

2(1 + kf)K{kO^ ^ _ £o(2Ac-g)' ^ ^ ^ k'^)K{k' f that is A;i^ = 0.985 (14) 

The interface of this new profile is relatively sharper (the width of the interface is now 

proportional to ^_^jK(k'>) ^ ^^^/z^ — ^7^^) compare to the size of the homogeneous 
domains which has now double, see Fig. (2). Again, this new stationary profile turns out to 
be linearly instable with respect to an "antiferro" perturbation of period AXq-h- 

Thus these families of profiles and instabilities enable to describe the one dimensional 
coarsening as a cascade of doubling process, leading from a pattern of wave length Xq-h 
composed of domains separated by interfaces to a single tanh(^^^^^a;) interface separating 
two semi infinite domains. Each of these successive intermediate profiles can be described 
by an element of the above family of sohton lattice ^(x, /c^,2" x Xc-h)- We thus have a 
family of segregation parameter {k^}, which are determined by the implicit relations 

2(1 + kf)Kik:,r = _ eo{2- X^c-nr ^ ^222n_ (^5) 
We have found numerically for the first of them 
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^= F = 0.6869795924 


A;^Ag = 0.400 V-^o 


A;^^ = 0.9851675587 


A;fA^= 0.496250 V-£o 


kl= 0.99997210165 


A:^A^= 0.499990V-£o 


A;| = 0.9999999999027 


A:|A^ = 0.49999846 V-£o 



(16) 



We see that {k'n} converges toward A;^ = 1 (single interface case) ; meanwhile the amplitude 
of the modulation k^Af^ goes toward ^|eo|/2, as can be seen in the second column of the 
table above. For large n, we can conclude from the implicit relation (15) that the ratio of 
the domain size to the interface width characterized by K{k^) behaves as 7r2'^~^. Each of 
the stationary profiles 



(17) 



2"7r ' 2"7r 

for which the interface width is proportional to -y==^^^ (which tends to -y=j in agreement 
with tanh(^^^^^a;)), is identically destroyed by the Langer "antiferro" instability. 



IV. AN ANSATZ FOR THE ID COARSENING PROCESS 

A. Non-symmetric soliton lattice Profile 

In order to describe one step of the coalescence process, i.e. the dynamics that start from 
^„(a;) and ends with the profile ^„+i(a;) , we will use another family of equilibrium profiles 
[19], solutions of (G-L) equation, which write: 

a{a, k) - kl^i5{a, k)Sn{4x^, k) 
1 — k/ y/aSn[4:X—j^, k) 

where a(a, k) = -2fcVa+i+fc^ ^ and f3ia, k) = ''-^=^ r • 

((l+A:2)2-12fc2+2(a+fe2/a)(l+fc2))^ ((l+A:2)2-12fc2+2(a+fc2/a)(l+fc2)) ^ 

It is still a periodic lattice of interfaces, but now, the mean value of the order parameter 
is non zero (non symmetric case) . It is controlled by the parameter a ^ 1 : if a goes infinity, 
we recover the previous family of periodic profiles. 



B. Ansatz for the continuous interpolation between two stationary states 



If we choose a to be equal to l + Zc' (where k''^ — 1 — A;^), we can then construct symmetric 
profiles using the sum of two non-symmetric ones. Indeed, using Gauss' transformation (or 
descending Landen transformation [20]), which relates the soliton lattice of spatial period 
2A (and of modulus k) to the soliton lattice of period A (and of modulus ii = j^), we have 

1 - ^^{i/j{k, X - I) + i/j{k, X + f )) - kSn{2x^, k) (19) 
l-^^(V'(fc,a;-A)+V(fc,a; + A)) = {I - k')Sn{{Ax + 2\)^, ^i) (20) 

where '0(/i;, x) — -0(1 + /c , /c, x). Thus, we then can show from equation (19) that 
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K{k) 



\/b — k'^, , A, , A,, 

1 (V^(fc, x--)+ iPik, X + -)) 



kK{k)Sn{2x^^,k). (21) 
A 



This is the solution of the G-L equation of period 2A. Moreover, if we use the fact that 



we can write 



(1 - k')K{k)Sni{Ax + 2X)^^, /i) = 2//X(^)>Sn((2x + A)^^^, //). 



A 



A 



(23) 



Then, using relation (20), the solution of the (G-L) equation of period A can be expressed 

as 



K{k) 



1 - {i^ik, x-X) + ^{k, X + A)) 



2iiK{ij)Sn{{2x + A)^^-^, /x). (24) 

A 



So, we see that both the initial state ^'*(x, A;^_]^, 2"~^Ac_//) and the final state 
'^*{x,k^^,2'^\c_H) of a step of the coalescence process can be describe, modulo a phase 
shift, by the same function : 



$(a;, k, 0) 



AK{k) 
A 



1 {ip{k,x 



(t)/2)\)+^{k,x 



0/2)A)) 



(25) 



with k — kf^ and A = 2'^Xc-h- Therefore we can describe the coalescence by a transformation 
at constant segregation parameter k, while the degree of freedom 0, associated with the 
relative phase between the two profiles, evolves in time from to 1 according to the C-H 
dynamics. 

This non-symmetric lattice of interfaces can be interpreted as a periodic sum of alternat- 
ing single interfaces (kinks and antikinks). In the same spirit as relation (9), if one forget in 
the infinite sum every two out of four interfaces, one gets : 



ip{x) ^ [tanh(7rs(a; — 4 * p)) — tanh(7rs(a; — p + 1))] . 



(26) 



Then (see Fig. (3)) adding ip{x + 2) to ■ip{x) enables to recover relation (9), while, after a 
translation, adding ijj{x + 1) and tjj{x) gives the soliton lattice of double period, because of 
the cancelation of half of the interfaces (annihilation of kinks and antikinks). 

If we look at the time evolution of the profile $(x, k, 0), starting from the region = 0, 
we can transform the (C-H) equation into a phase field equation, replacing ^ $(a;, k, 0) by 
^$(a;, k, 0(t)).^. The dynamics will be similar to a spinodal decomposition, with growing 

and saturating exponentially. ^$(a;, k, 0) is the most unstable mode founded in Langer's 
linear stability analysis, characterized by the alternate growth and decrease of domains 
("antiferro" mode). Note that when Langer was studying the most instable perturbation, 
he was looking at the hnearized version of C-H equation around ^*{k, x) — ^{x, k, Xc-h) ■ 



= + 6^*2 _ <^ = + n X (n + 1)**' - V') if. 



(27) 
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Cip = Eip is the Lame equation, for n = 2 (here Eq = 1). This equation doesn't have 
simple (algebraic) exact eigenfunction of period 2Xc-h [21]- ^ ^{x,k,(f)) for = is not 
an exact eigenfunction either [22]. Nevertheless, it happens to be a good approximation 
for the eigenfunction of lowest eigenvalue. Due to the concavity of J^(0) around = (see 
below Figure (5)), this eingenvalue will be negative, triggering a linear destabilization and an 
exponential amplification of the perturbation, i.e. an exponential growth of the translation 
(j) with time. 

hanger's phenomenon of "antiferro" instability appears due to the existence of two pos- 
sible directions for displacement of the interfaces " tanh " (or of the non-symmetric lattice of 
interfaces ■0), one with a positive velocity (+^) and one with a negative one (— The four 
different kinds of interfaces present in a cell of length 2Xc-h have alternately a positive or 
a negative velocity. This can be seen as the existence of two antisymmetric patterns [23] , or 
building blocks for the leading instability around a intermediate state ^*(a:, k^, 2"^^ x Xc-h) 
(see Figure (4)). ±-^ip{x), these two building blocks, are associated with the two pairs of 
interfaces, ip{x) and ip{x + 2) which have been used to construct our ansatz. 

Note that in Langer's analysis, the breaking of symmetry for the choice of the antiferro 
cell, corresponds here to the freedom we have when choosing the range of variation of 6 : we 
could have chosen to go from to —1, ending after a step of coarsening with the symmetric 
pattern, or equivalently, a pattern translated of half a period. 

C. Energy Landscape 

In order to prove the usefulness of this ansatz, we have plot the energy averaged over the 
final period, J-'{4>) = J F{^{x, k, (f)))dx, as a function of the parameter 0, keeping k constant. 
We see for example in Figure (5) that the value = correspond to a local maximum of 
energy, while 0=1 (or —1) is a minimum. Note that there is no energy barrier in this 
particular energy landscape, in agreement with linear stability analysis. 

D. Approximation for adiabatic evolution with constant k 

Relation (15), enables to write an implicit relation between kn and : 

^l + k^^^,K{K+i) = 2./l^K{kl) (28) 

It is in fact different from the relation (22) obtained by the Landen transformation. Never- 
theless in the region k > k^ — 0.687, the two transformation almost coincide, as can be seen 
on Figure (6). As this is especially true close to k = 1, and as this region is rapidly reach 
after the second or third iteration, we pretend that the process of coarsening can be describe 
with a reasonable accuracy by our antsatz at constant k. A slight change of segregation 
parameters during the n*'* doubling process, from k = Landen{k^) when = to /c = k^_^^ 
when = 1, is present in the dynamics, as seen in the following table 



A:^ = 0.6869795924 


Landen-^{k^ = Gaussikl) = 0.7070743852 


L(ui(h:i){l4) = 0.9826316738 


A-^ = 0.9851675587 


Landen{kf) = 0.9999720868 


^= 0.99997210165 


Landen{k'^ = 0.9999999998 


A;| = 0.999999999902745 



(29) 
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But it has only a minor effect in terms of profile shape /c,0). This means that the 
dynamics of the parameter k, which affect the value of the modulus by 0.25% during the 
first step of coarsening, becomes negligeable as k goes closer to 1. So we can conclude that 
this " /c-dynamics" is irrelevant ; this parameter can be considered as constant during the 
evolution of (p. 

V. DISCUSSION ON THE HYPOTHESIS 

Our analytic method rely on the assumption that at each step of the dynamics, the 
system can be characterized by a specific spatial period : we need therefore to discuss 
how this approach is relevant to the general case where noise is present. We have noted 
numerically that, for the spinodal decomposition, the average size of the modulation is Xc-h, 
with a deviation of less than one percent from the value predicted by the linear theory. It 
does not mean that, in a real system, each domain has a length scale of \c-h, but that the 
distribution of the domains' length will be centered around Xc-h [6]. The coalescence events, 
due to initial fluctuations in the periodicity of the pattern selected just after the quench, 
can be neglected during the initial growth of the amplitude of the modulation. Only after 
this initial growth has saturated (as can be seen from Figure (1)), starts the coalescence 
process, which will then dominate the dynamics : the typical length scale of the structures 
increases slowly with time. 

As suggested by [1,15], we can suppose that during the ideal coalescence process, each 
lattice of interfaces will experience an antiferro instability. By ideal coalescence, we mean 
a process which breaks as few symmetries as possible. But in a real system, this instability 
will concern a region of finite size, where it choose a certain sublattice, or a range for (f) (for 
example, varies from to 1), while it is the opposite choice in the neighboring region ( (j) 
varies from to -1). Thus on the overall, the global symmetry is recover. During each step 
of the process, the width of the domains will locally double. But due to non synchronization 
between regions, for the system as a whole, the average length scale will vary continuously. 

These results indebted to Christophe Josserand for many fruitful discussions and com- 
ments together with numerical help. 
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FIGURES 



FIG. 1. Time evolution of the order parameter ^{x,t) for e = —1, dx = 0.1227. (a) initial 

conditions at t = are taken randomly with a very low amplitude (5 • 10~^) ; (b) at time t = 15, 
the amplitude of the small scale spatial modulations has been damped by the (C-H) dynamics, 
while only long wavelength contributions are still present. ; (c) at t = 225, the spatial modulation 
has almost reached its final amplitude, keeping roughly the same number of peaks as before ; (d) 
at t = 1800, we observe that the number of domains has decreased from the coarsening dynamics. 

FIG. 2. Profiles of the two first metastable solutions of the (C-H) dynamics, with kl = 0.687 
and A;2 = 0.985, corresponding to the first coarsening process. 

FIG. 3. Construction of the two first steady solutions of the (C-H) dynamics, with = 0.687 
and /c| = 0.985, using a superposition of the non-symmetric profile ipik, x), itself stationary solution 
of the (C-H) equation. By changing the phase shift between the two profiles entering into the linear 
combination, one obtains two different symmetric profiles, of periods A and 2A. 

FIG. 4. Langer's most instable perturbation mode of destabilization of the soliton lattice is 
identified with ^ ^{x,k,(f)) at ^ = 0. It is composed of two antisymmetric patterns, plotted 

in dotted (plain) line, evolving toward right (left) at velocity +^ ("^)' causing an "antiferro" 
instability leading to a period doubling of the pattern. They are the spatial derivative of the initial 
non symmetric profile ip{x) which has been used to construct our ansatz in Figure 3. 

FIG. 5. Profile of the free energy landscape during a coarsening process, F(^). It starts at 
= for a configuration characterized by the segregation ratio A;f = 0.687 for which the energy per 
unit length is F(^ = 0)~ —0.135; one sees that in this region, the free energy is a concave function 

of (f) and thus, the associated pattern is linearly instable. The elementary step of the coarsening 
process ends for (f) = \ associated with a pattern characterized by the segregation ratio ^2 = 0.985 
for which the energy per unit length is F(^ = 1)~ —0.45. In the region (j) = 1, the free energy is a 
convex function of (f). 
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FIG. 6. Comparison between the Landen transformation (upper solid line) and 
the implicit relation between consecutive stationary steady states of (C-H) equation 
(1 + k'^_^i)2 K (kn+i)=2{l + k'^)2K(kn) (lower dash line) in the region between k = kf = 0.68 
and A; = 1, corresponding to the region of interest for the coarsening process. The Landen trans- 
formation relates the segregation parameter of a soliton lattice of period A with the segregation 
parameter associated with a lattice of period 2A. It is the generalization for the cnoidal function of 
the usual relation sin 20 = 2 sin cos 0. The circles mark (from left to right) /cq, Landen""^ (fef), /cf , 
and LanderT^ {k2) ■ If the doubling process associated with the coarsening were only a phase shift 
of between and 1, the two curves would have coincided, i.e. the implicit relation (28) would 
be equivalent to the Landen transformation. As it is not the case, there is a A; change during the 
doubling process, but only of a few percent, as can be seen on the figure above. Moreover, one sees 
that, as the coarsening process takes place, k reaches values closer and closer to 1 for which the 
disagreement becomes negligible : the k change becomes smaller and smaller. 
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